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Abstract 

We have developed a new technique for inclusive reconstruction of the energy 
of B hadrons. The excellent efficiency and resolution of this technique allow us 
to make the most precise determination of the 6-quark fragmentation function, 
using e+e^ — > decays recorded in the SLD experiment at SLAC. We 
compared our measurement with the predictions of a number of fragmentation 
models. We excluded several of these models and measured the average scaled 
energy of weakly-decaying B hadrons to be {xb) = 0.714 ± 0.005 (stat) it 
0.007 (syst) ± 0.002 (modebdependence). 
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In high-energy strong-interaction processes quarks and gluons are not observed directly, 
but appear as jets of colorless hadrons. This fundamental process of 'jet fragmentation' 
affects all high-energy physics measurements involving strongly-interacting particles, but 
is only poorly understood at a quantitative level. The fragmentation of heavy quarks is of 
particular experimental interest since many expected signatures of new heavy particles, such 
as Higgs and SUSY particles, involve decays to b quarks. It is hence vital to understand the 
production and properties of 6 jets. Here we present a significantly improved determination of 
the 6-quark fragmentation function, D{xb) = {l/o')da/dxB, measured using bb decays, 

where xb = Eb/ Ebeam represents the fraction of the 6-quark energy retained by the weakly- 
decaying B hadron. 

In Quantum Chromodynamics (QCD) the 6-quark mass serves as a cutoff for collinear 
gluon radiation. The distribution of 6-quark energies prior to hadronisation can there- 
fore be calculated perturbatively |[l|-|5|. However, the additional effects which yield the 
experimentally-accessible distribution D{xb) are non-perturbative, and have been studied 
phenomenologically via a number of different approaches [|T]-0]. Measurements of D{xb) 
provide direct tests of these perturbative QCD and model predictions. 

D{xb) has been measured previously by reconstructing the energies of B hadrons 

that decay semi-leptonically {B IvDX) in — > hb events. In these studies the inclusive 
B selection efficiency was smaller than 1%, and much lower for Eb < 20 GeV. The resulting 
low-statistics samples, and hmited energy resolutions, yielded poor constraints on the shape 
of the distribution. 

We have developed a new technique for measuring Eb using only kinematic information 
from charged tracks. Our 307 Mpixel CCD-based vertex detector, combined with the micron- 
sized SLC interaction point (IP), allows us to reconstruct 5-decay vertices and the B- 
hadron flight direction very accurately. The method yields a significantly higher 5-selection 
efficiency and superior energy resolution, both of which are almost independent of Eb-, and 
has very low bias since it does not use a beam-energy constraint. This allows us to measure 
the shape of D{xb) with sufficient precision to make stringent tests of 6- fragmentation 
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model predictions, reduce model-dependent systematic errors, and therefore discriminate 
among these models for the first time. Furthermore, our technique for reconstructing B- 
hadron energies is of direct relevance to studies of other important properties of heavy-quark 
systems, such as i?-hadron lifetimes and neutral i?-meson mixing. 

We used 150,000 hadronic decays produced in e~^e~ annihilations at the SLAG Linear 
Collider (SLC) and collected in the SLC Large Detector (SLD) between 1996 and 1997. A 



description of the SLD can be found elsewhere |117| , [18| . The trigger and selection criteria 



for hadrons events are described elsewhere |jT6|. This analysis used charged tracks 



measured in the Central Drift Chamber (CDC) ||T9| and in the upgraded CCD Vertex De- 
tector (VXD) pO|, with a momentum resolution of crpj^/p± = 0.01 © 0.0026p^, where p± is 
the track transverse momentum with respect to the beamline, in GeV/c. The centroid of 
the SLC IP was reconstructed with a precision of 4.4/im (30/im) in the plane transverse to 
(containing) the beamline. Tracks from identified 7 conversions and or A° decays were 



removed from consideration, and only well-reconstructed tracks [|T^ were used for i?-hadron 
tagging and energy reconstruction. 

Weakly-decaying B hadrons were identified by exploiting their long lifetimes and large 
masses relative to light-flavor hadrons. Each hadronic event was divided into two hemi- 
spheres by the plane containing the IP and the normal to the thrust axis. A topological 



vertexing algorithm pl[ was optimized for this analysis and applied to the set of well- 
reconstructed tracks in each hemisphere in an attempt to reconstruct a i?-decay vertex. 
Vertices were required to be separated from the IP by at least 1 mm and to contain at least 



two tracks. A candidate vertex was found in 32,492 hemispheres p2 . 

In each hemisphere the total energy, Ech, momentum, Pch, and invariant mass, Mch, of 
the vertex-associated tracks were calculated by assigning each track the charged-pion mass. 
Due to neutral decay products and tracks missed from the vertex, P^h can be acollinear with 
the true B flight direction, which we estimated independently of the track momenta by the 
unit vector v along the line joining the IP to the reconstructed vertex position. Our accurate 
knowledge of the IP and vertex positions yielded an excellent angular resolution of 26 mrad 



on V. This allowed us to calculate the net transverse momentum, relative to v, of the particles 
from the B decay that were not associated with the vertex, Pf = Pch — {Pch'V) v, and hence to 



improve our estimate of the mass of the decaying particle via Mp^ = y M^f^ + Pf + \Pt\] this 
quantity is a strong discriminator for selecting B hadrons [E^. We required 2.0 GeV/c^ < 



< 2 X Mch to select 19,604 candidates, estimated to be 98.2% B hadrons, the main 
background (1.6%) being charmed hadrons in cc events. The efficiency for selecting a 

true B hadron is 40.1%. For this sample, on average 92% of the reconstructed true S-decay 
tracks were associated with the vertex, and 98% of the vertex- associated tracks were from 
true B decays. 



The reconstructed _B-hadron energy, E^'^, can be expressed as \j Mq + P^ + P^i + Ech^ 
where the combined mass of the missing particles, Mq, and the missing momentum along iT, 
Pqz, are the only unmeasured quantities. If we assume a 5-hadron rest mass, Mb-, we can 
eliminate one of the two unknowns and calculate an upper bound on Mq: 



M^ < Mi - 2Mb ^M2, + P? + Ml ^ Mo^^, , (1) 

where equality holds when Pqi = in the S-hadron rest frame. Since, from phase-space 
requirements, small values of Pqz are the most probable, and the average -B-decay multiplicity 
is high, the true Mq tends to be close to M^max Hll. We assumed Mb = 5.28 GeV/c^, equated 



Mq with the measured Mo^nax, and calculated E^'^. This estimate of the B energy is best 
when the -B-decay kinematics are well constrained by the vertex-associated tracks, i.e. when 
MQmax — p5|. Also, the small non-i? background is concentrated at large Momax- In 



order to improve the energy resolution, and reduce further the background, we required 
—1 < Mq^„^ < [1.1 + f{E^'^)]'^, where the explicit i^^^'^'^-dependence was chosen so that 
the B selection efficiency is only weakly energy-dependent. The efficiency is above 3% for 
Eb > 10 GeV; the average value is 3.9%. 1920 candidates were selected, with an estimated 
B purity of 99.3%. 

We examined the distribution of the normalized difference between the reconstructed 



and true 5-hadron energies, {E^b'^ — E^b^^) / E^b^^ . This resolution can be characterised ||24 



by a double Gaussian function centered at zero. We found that the narrower Gaussian 
represents 83% of the population and has a width of 10.4%. This resolution depends only 
weakly on E'j^^'^; in particular it remains better than 15% even for B energies close to the 
mass threshold, which is a significant advantage of this technique. The estimated 0.7% non- 
B background was subtracted bin-by-bin from the distribution of the reconstructed scaled 
i?-hadron energy, D^'^'^{x^^'^), which is shown in Fig. 1. 

We tested several models of 6-quark fragmentation. These models are formulated in 
terms of experimentally inaccessible variables and must hence be implemented in an iter- 
ative fragmentation algorithm in order to derive the measurable quantity D^'^^{x^j^^). We 
employed our JETSET 7.4 [^-based Monte Carlo simulation program [lT6| , p3| to gener- 



ate e+e bb events according to each model considered. 5-hadron energies were re- 
constructed, according to our algorithm, from the fully-simulated event sample to derive 
which was compared with the data using a binned x^- The was minimised by 
repeating this procedure under variation of the input parameter(s) of the model. The fitted 
model predictions and best values are shown in Fig. 1. 

Within this context the Bowler Kartvelishvili et al. |^ and Lund [^] models 

reproduce the data; the models of Braaten et al. M, Collins and Spiller ||T2[ and Peterson 



et al. ||TOl have a confidence level less than 0.1% and are not consistent with the data. 



We also tested the Monte Carlo models HERWIG 5.7 [||] and UCLA [H], which contain 
no explicit free parameters to control D{xb)- The UCLA model is consistent with the data 
and the HERWIG model is not (Fig. 1). 

In order to allow other models to be compared with our data, e.g. those of |l|-0,||-0, we 
corrected for the effects of the selection and reconstruction procedures that were applied. 
We estimated the true weakly-decaying i?-hadron scaled-energy distribution, D{xb), from 
the (background-subtracted) reconstructed distribution, D^^'^{x'^^^); for each bin i: D{xB)i = 
EfcMjfcD'"^'^(x^'^)fc/ej. The selection efficiency e and the unfolding matrix M were calculated 
from the simulation using in turn each of the four fragmentation models that were consistent 
with the data, with the respective fitted parameter values. This unfolding procedure assumes 



a smooth underlying distribution and is explicitly model dependent, but we quantify this 
(below) using the variation among the four resulting distribution shapes. 

This set of consistent models is small and does not appear to span the range of potentially 
acceptable shapes. We therefore considered a number of ad hoc functional forms for the true 
xb distribution, and found four that yielded a reconstructed distribution consistent with 
the data in Fig. 1: an 8th-order polynomial, the Peterson function, and two generalisations 



of the Peterson function |jT^. The respective parameter value(s) was (were) optimised in a 
manner similar to that described above. Each of these four fitted functions was also used to 
calculate e and M from the simulation and to derive D{xb)- 

Each of the eight unfolded distributions was normalized to unit area, and in Fig. 2 we 
show their binwise average; the band includes the r.m.s. deviation. This represents our 
best estimate of the true B energy distribution D{xb)- The corrected distribution is, by 
construction, smoother than the measured distribution, and the band provides an envelope 
within which acceptable predictions should fall. The constraint on the shape is much stronger 
than in any previous determination |T^-|TB|. 



An important advantage of our method is that it is not biased by tracks that were 
not reconstructed or attached to the vertex P3. We considered other potential sources of 



systematic uncertainty which derive from the modelling of our detector response. In each 
case the simulated events were reweighted or adjusted in order to reproduce the variation 
in question, and the entire analysis was repeated. Variation of quantities such as the track 
reconstruction efficiency, the point resolution of the vertex detector and the momentum and 



dip angle resolutions of the tracking system within their uncertainties affects the B 
selection efficiency and, in some cases, the energy resolution. However, there is little effect 
on the shape of the distribution or on the mean value, (xb)- In no case was the of 
a model or function test changed significantly. The largest effect on (xb), ±0.005, arose 
from the uncertainty in the momentum resolution, which was measured in the data using 
e'^e" — i> fi^fi~ events. An ad hoc correction was applied to the simulation to reproduce the 
measurement, and the full effect of this correction was considered as a symmetric systematic 
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uncertainty. As a cross check we changed the upper cut on Mq^^^ to fixed values between 
1 and 5 (GeV/c^)^; the change in {xb) was smaller than the statistical error. 

We also considered the uncertainties on a large number of measured quantities related 
to the production and decay of B and charmed hadrons. These are used in the simulation 
and potentially affect the values of e and M used to unfold the data. We varied each quantity 
by the error on its measured value; none of these variations affects our conclusions. The 
production of primary excited B hadrons, collectively denoted B* and -B**, which decay 
into the weakly-decaying B hadrons that we measure, was investigated in more detail, and 
we varied independently the fraction of primary B* (B**) from zero to unity (0.5). No 
significant effect on the shape of D{xb), or on the values for the model comparisons, 
was observed. The largest effect on (xb), of ±0.002, was due to variation of the number of 
charged tracks and K° produced per S-hadron decay. In each Xb bin the sum in quadrature 
of the detector- and physics-related systematic uncertainties on D{xb) is much smaller than 
the statistical error; they are not shown in Fig. 2. 

From the eight distribution shapes that are consistent with our data we extract the mean 
value of the scaled energy of weakly-decaying B hadrons in Z° decays: 

(xb) = 0.714 ± 0.005(stat.) ± 0.007(syst.) ± 0.002(model). 

This is the most precise of the world's measurements that take the D{xb) shape dependence 
into account, and the uncertainty is relatively small since we have excluded a wide range 
of shapes. Our result is consistent with a recent average over many measurements 
of (xb) = 0.702 ± 0.008 ||3^. We also calculated the second and third moments to be 
{xD = 0.537 ± 0.011 (exp.) ±0.003 (model) and = 0.417 ± 0.012 (exp.) ±0.004 

(model). 

In order to derive results for the inclusive sample of primary B hadrons, one must assume 
values for the production fractions of B* and B** mesons. Postulating a B* production 
fraction of 0.75, expected from naive spin counting, leads to (^^ bY"™'"''^^ = 0.718. Postulating 
independently a B** production fraction of 0.25 yields {xbY""^""^^ = 0.728. 



In summary, we have developed a new, inclusive technique for reconstructing the ener- 
gies of B hadrons. It has substantially higher efficiency and better energy resolution than 
previous methods. We have employed this technique to measure the scaled-energy distri- 
bution of weakly-decaying B hadrons produced in e'^e" — > decays with unprecedented 
precision over the entire kinematic range from the B mass to the beam energy. As a result 
we are able to exclude several models of 6-quark fragmentation, including the widely-used 
JETSET+Peterson model. 

We thank the personnel of the SLAC accelerator department and the technical staffs 
of our collaborating institutions for their outstanding efforts on our behalf. We thank 
M. Bowler and M. Seymour for helpful comments. This work was supported by the U.S. 
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FIGURES 
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FIG. 1. Distribution of the reconstructed scaled energy of weakly decaying i?-hadrons (points); 
the errors are statistical. The predictions of eight models are shown as histograms. 
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FIG. 2. Unfolded distribution of weakly-decaying scaled 5-hadron energy (points). The errors 
are statistical only and do not include point-to-point correlations. The band represents the envelope 
of acceptable functions (see text). Also shown (circles) is the best previous measurement [15]. 
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